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Abstract. We present the Geodesic Deviation Equation (GDE) for the Friedmann- Robertson 
Walker (FRW) universe and we compare it with the equation for Bianchi type I model. We 
justify consider this cosmological model due to the recent importance the Bianchi Models have 
as alternative models in cosmology. The main property of these models, solutions of Einstein 
Field Equations (EFE) is that they are homogeneous as the FRW model but they are not 
isotropic. We can see this because they have a non-null Weyl tensor, which is zero for FRW 
model. We study some consequences of this Weyl tensor in the GDE. 



1. Introduction 

Bianchi cosmological models are homogeneous solutions of EFE. These models generalize FRW 
model, although they are homogeneous, as FRW one, generally they are anisotropic. Standard 
Cosmological Model has survived almost all observational tests [8J, however Bianchi cosmologies 
haven't been discarded. Today the Universe is highly isotropic, but in early times it couldn't 
been so, and it is an open question how anisotropics died. Standard Model fails to explain 
quadrupole anomaly [5j, obtained from Cosmic Microwave Background data. This anomaly 
have a plausible explanation with Bianchi models, for that reason these models are considered 
today [5J. Thus we considered the idea of studying Geodesic Deviation Equation(GDE) in some 
of these models. 

2. 1+3 Covariant description 

We assume a fluid description for the cosmological model. The 4- velocity u a of preferred world- 
lines is normalized, u a u a = — 1. Given this 4- velocity, we define the projection tensor h a p, the 
expansion scalar O, the shear cr a p and the vorticity [7], [2]. In FRW case cr a p = uj a p = 0. 
These quantities are related with Weyl tensor. In FRW case Weyl tensor is zero, so it is confor- 
mally flat, but it is not the case in Bianchi models. In the 1 + 3 covariant description we can 
split the Weyl tensor into an electric E a @ and a magnetic part H a p [2] , [7] . 



3. Deviation Equation 

In a space-time torsion free, if V a is the normalized tangent vector field of a fiducial geodesic, 
parameterized by an afflne parameter is, with V a V a = £, with £ = +1,0,— 1 if the geodesies are 
spacelike, null or timelike respectively; if r] a is the deviation vector for the congruence, then the 
evolution of the separation or the deviation vector r] a evolves according to GDE [lj: 
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where ^" — V^ 7 V 7 T 6 a " is the covariant derivative along the geodesic. 

Now, let's see how is this equation in Bianchi type I model, given by the metric[3], 

ds 2 = dt 2 - a 2 {t)dx 2 - b 2 {t)dy 2 - c 2 {t)dz 2 (2) 

If we define r{t) := a(t)b(t)c(t), from the energy-momentum tensor conservation we get an 
analogous expression to a known Friedmann equation [3] , [I] : 
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where is the Energy density and p is the pressure. As we see in the GDE, there is a dif- 
ference if space-time has a nonzero Weyl tensor. Ellis and van Elst in [lj consider this equation 
for FRW space-time. 

For Bianchi type I model vorticity and Magnetic part of Weyl tensor are equal to zero, 
— 0, H a p = 0. The Electric and shear tensors are non-null and diagonal. For that reason, 
if we define F = E^V^rf 1 and E — —V a u a , if we suppose a perfect fluid, q a = 0, ir a p = we 
have that GDE for this space is: 
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where A is the cosmological constant. The first part of this expression is the force term 
representing the perfect part of the fluid [2J and only gives isotropic deviation. In FRW case we 
only have this term and reflects spatial isotropy of space-time, deviation vector only changes in 
magnitude but not in direction. The other terms change the direction of r] a . 

If we suppose an equation of state p = we, (w = for dust, w = | for radiation), then, for 
A = weget[3]: 
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Given r as a function of t we can get a^t), b(t) and c(t). The electric Weyl tensor tends 
to zero. It is diagonal and its components decay. For this model the vorticity and the Mag- 
netic part of Weyl tensor are null. The EFE imply that q a = 0. This means we have a fluid 
without viscosity. However, we can have an imperfect fluid because the model allows to have an 
anisotropic pressure 7r a p. But this anisotropic pressure tensor must be diagonal, following EFE's. 
Given r as a function of t we have expressions for a(t), b(t) and c{t). The evolution of the three 
scale factors for w = is shown in Figure 1 and for w = | in figure 2. We can see that in radi- 
ation case anisotropy is greater than in matter dominated case, but in both cases anisotropy dies. 

The a scalar defined by a 2 = \cr a ^a a p also measure the degree of anisotropy in the model. 
In Figure 4 we illustrate the evolution of a 2 . 




4. Conclusions 

The evolution of the Electric Weyl tensor show us that this model tends to isotropy as t grows, 
the model isotropizes. Also the evolution of shear scalar. In FRW only the magnitude 77 will 
change along a geodesic, while its spatial orientation remain fixed [lj. So, in early times the 
deviation is anisotropic, the orientation is important, for example, along the z axis deviation will 
grow, while along x and y axis the effect of Weyl Tensor is to contract the deviation, because 
there are the other terms that are also present in FRW. Here we considered a perfect fluid, 
^a(3 — 0. Along x axis, the effect of anisotropy is minimal. 
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Figure 4. Electric Weyl Figure 5. Shear scalar 

tensor for A = |, w = 5. evolution for A = |. 
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